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Abstract- A  broadband  time  domain  reverberation  and 
clutter  model  has  been  developed  to  aid  in  the  development  of 
a  better  understanding  of  the  environmental  and  system 
characteristics  which  drive  sonar  clutter  in  inhomogeneous 
waveguides.  The  theoretical  background  of  the  model  is 
presented  and  its  ability  to  model  both  diffuse  background 
reverberation  levels  and  clutter  characteristics  for  a  site  on  the 
Malta  Plateau  south  of  Sicily  are  evaluated. 

I.  Introduction 

Sonar  clutter  in  littoral  environments  increases  the 
probability  of  false  alarm  for  a  given  probability  of 
detection.  This  is  because  clutter  adds  to  the  length  of  the 
tails  of  the  reverberation  envelope  probability  distribution 
function  (pdf),  moving  the  statistics  away  from  the  Rayleigh 
canonical  form.  Clutter  can  be  caused  by  target-like  features, 
either  natural  or  man-made,  or  by  non-Gaussian 
distributions  of  the  scatterers.  Typically,  high  bandwidth 
and/or  highly  directive  systems  have  more  problems  with 
clutter,  because  as  the  size  of  the  scattering  patch  is  reduced, 
the  pdf  of  the  generally  non-Gaussian  scatterer  distributions 
becomes  resolved  by  the  system.  Recent  work  [1]  has 
demonstrated  that  for  single  path  saturated  propagation  to  a 
collection  of  exponentially  sized  scatterers,  the  scattering 
envelope  has  the  analytic  form  of  the  K  distribution,  in  the 
physical  optics  limit.  The  K  distribution  shape  parameter, 
which  controls  the  degree  of  non-Rayleighness  of  the 
scattering  envelope,  was  shown  in  turn  to  be  directly 
proportional  to  the  number  of  scatterers  within  the  sonar 
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footprint.  For  scatterer  numbers  greater  than  approximately 
40,  the  scattered  envelope  was  shown  to  revert  to  Rayleigh. 

For  reverberation  in  inhomogeneous,  range  dependent 
shallow  water  waveguides,  it  becomes  necessary  to 
generalize  the  results  of  Abraham  and  Lyons  to  obtain 
predictions  of  interest  for  notional  sonar  systems.  In  range 
independent  waveguides  the  propagation  physics  along  with 
the  system’s  directivity  and  bandwidth  control  the  number 
of  scatterers  that  contribute  to  the  reverberation  envelope 
statistics.  For  example,  for  rough  surfaces,  the  shallow  angle 
paths  which  persist  the  longest  due  to  low  attenuation,  and 
travel  the  fastest  due  to  the  higher  group  speeds,  generally 
also  generate  the  least  amount  of  scattering  locally  from  a 
given  roughness  distribution.  Conversely,  the  heavily 
attenuated  steeper,  slower  paths  generate  relatively  more 
scattering  locally  for  a  given  level  of  incident  intensity. 
This  means  that  there  is  a  trade-off  of  the  different  paths  as 
a  function  of  range,  with  the  steeper  paths  being  dominant  at 
short  range  and  the  shallower  paths  becoming  more 
important  at  long  range  [2].  The  dispersion  of  the 
waveguide,  which  leads  to  time  spread,  together  with  the 
bandwidth  determine  the  degree  to  which  the  various 
multipath  contributions  to  the  reverberation  are  independent. 
The  number  of  range-wise  scatterers  contributing  to 
reverberation  is  controlled  by  the  bottom  critical  angle,  the 
vertical  directivity  of  the  source,  the  scattering  angle 
dependence  of  the  scattering  kernel,  the  loss  vs.  angle 
relation  of  the  bottom  interaction  loss,  and  the  bandwidth, 
while  it  is  the  horizontal  beamwidth  alone  that  controls  the 
number  of  independent  azimuthal  scatterers  contributing  to 
the  reverberation  for  the  monostatic  systems  considered 
here. 
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In  order  to  incorporate  the  above  considerations  into 
developing  a  better  understanding  of  the  environmental  and 
system  drivers  of  sonar  clutter  in  inhomogeneous 
waveguides,  a  broadband  time  domain  monostatic 
reverberation/clutter  model  has  been  developed  based  on  the 
R-SNAP  model  [2].  The  R-SNAP  model  is  based  on  the 
theory  of  coupled  normal  modes  and  exploits  the  knowledge 
of  modal  group  speeds  and  dispersion  to  generate 
approximate  reverberation  time  series  using  the  narrow 
band  approximation.  Until  now  the  model  has  been  used 
solely  to  predict  the  reverberation  intensity  of  monostatic 
systems  as  a  function  of  the  environmental  and  system 
parameters,  although  it  was  also  capable  of  predicting  the 
backscattered  levels  returned  from  deterministic 
environmental  discontinuities  in  the  submarine  environment. 
Here  the  model  has  been  extended  to  generate  reverberation 
time  series  consistent  with  the  closed  form  intensity 
estimates,  generated  in  a  Monte-Carlo  manner  from 
realizations  of  the  bottom  scatterers.  In  this  paper  we 
describe  the  theoretical  background  of  the  model  and 
evaluate  its  ability  to  model  both  diffuse  background 
reverberation  levels  and  clutter  characteristics  for  a  site  on 
the  Malta  Plateau  south  of  Sicily. 


II.  Theory 


A.  Reverberation  Simulation 

The  R-SNAP  model  is  based  on  the  C-SNAP  coupled 
normal  mode  propagation  model  [3].  In  C-SNAP  the  total 
field  at  a  scatterer  depth  zscat  is  given  as  a  function  of  the 
range  dependent  complex  modal  amplitude  an  and  the 
range  dependent  mode  depth  functions  (j)n 


P{r’  2  scat  )  =  Jla„(r)t(/,Zscat)> 


n= 1 


(i) 


In  R-SNAP  we  are  instead  interested  in  the  field  incident 
on  the  scattering  plane  at  the  scatterer  depth.  This  is 
accomplished  by  decomposing  the  mode  shape  functions  (/) 
into  up  and  down-propagating  plane  waves  c/)~  [4].  The 
field  incident  on  the  scattering  plane  at  depth  zscat  is  then 


~{PZscat)  = 


N 

n= 1 


): 


(2) 


and  the  field  backscattered  from  an  array  of  objects  at 
ranges  rf  with  a  scattering  function  ss[r^On^Om)  is  then 
given  as 


-i  ^ss(r\Gn,G^(l(r\zJjam(r)drr 


(3) 


Equation  (3)  is  an  expression  for  the  scattered  field  in  the 
frequency  domain.  Since  reverberation  is  observed  in  the 
time  domain,  the  Fourier  transform  of  (3)  is  required  to 
perform  a  simulation.  This  Fourier  transform  can  be  well 
approximated  in  the  case  of  relatively  narrowband 


excitation  through  the  use  of  the  narrowband  approximation 
for  coupled  propagation  in  range  dependent  media  [5].  For 
one-way  propagation  the  narrowband  approximation  for 
range  dependent  waveguides  is 
P(t,r,zr) 


=  2Re^ 


i 


V1/A °>  -dn(r)r 


exp 


-icot  - 


2(l/ Act)2  -dn  (r)r) 


(4) 


where  CO  is  the  frequency  at  which  the  solution  to  (1)  was 
obtained,  Sn  is  the  first  derivative  of  the  mode  n 
wavenumber  kn  with  respect  to  circular  frequency,  dn  is 
the  corresponding  second  derivative,  and  A  CO  is  the 
bandwidth  of  the  source  signal  in  rad/s. 


B.  Scattering  Strength  Decomposition 


In  the  work  [1]  on  single  path  scattering  statistics,  the 
scattering  strength  of  the  spheres  with  exponentially 
distributed  radii  was 

ss(r)  =  al  2, 

(5) 

the  geometric  optics  limit  for  a  rigid  scatterer.  For  rough 
surface  scattering,  we  argue  that  the  scattering  strength  may 
be  decomposed  into  the  product  of  an  angularly  dependent 
term  and  a  spatial  distribution.  Whereas  in  [1]  the  spatial 
distribution  was  discrete  with  sufficient  spatial  randomness 
to  guarantee  phase  saturation,  in  rough  surface  scattering  we 
are  dealing  with  a  continuous  distribution.  Thus  our  model 
is 

ss(r,6n,6a)  =  Ti{r)f(a),6n,6m,'E), 

(6) 

where  T]  is  the  height  profile  of  the  roughness,  6 n  m  are  the 
incident  and  scattered  grazing  angles  associated  with  modes 
n  and  m ,  and  E  is  a  vector  of  environmental  information 
containing  information  about  the  water  and  bottom  sound 
speeds,  densities  and  attenuations.  In  this  case  it  may  be 
shown  that  under  the  first  order  perturbation  approximation 
[6] 


Ano.A.  e) 


where  B  is  the  boundary  condition  operator  for  the 
unperturbed  boundary  described  by  E  at  the  rough  surface, 
b  is  an  associated  rotational  operator,  and  p{kn^j  is  a 
vector  containing  a  unity  amplitude  plane  wave  incident  on 
the  mean  scattering  surface  at  the  grazing  angle  associated 
with  mode  n ,  along  with  the  associated  reflected  and 


transmitted  components  for  the  unperturbed  boundary 
described  by  E  . 


An  analogous  description  for  arbitrary  scattering  strength 
angular  dependence  may  also  be  utilized.  As  long  as  the 
correlation  length  scale  of  the  scatterer  distribution  is 
smaller  than  half  an  incident  wavelength,  a  very  good 
approximation  is 

f{wOnA>  E)=l(f(^)/20, 

(8) 

where  SS  (  CO,  6n ,  6m  )  is  some  parametric  form  of  the 
scattering  strength,  derived  from  data  for  instance. 

III.  Malta  Plateau  Data 

A.  Reverberation  Characteristics 

The  Malta  Plateau  is  a  shallow  (70-200  m)  site  in  the 
Straits  of  Sicily  north  of  Malta.  It  is  characterized  by  a 
muddy  bottom  lying  over  consolidated  sediment/rock  which 
becomes  exposed  in  a  ridge  running  north-south  between 
Malta  and  eastern  Sicily.  In  Fig.  1  the  sediment  thickness  to 
the  consolidated  sediment  is  illustrated  in  color  with  the 
bathymetric  contours  superimposed. 


Figure  1.  Sediment  thickness  and  bathymetry  on  the  Malta  Plateau. 

In  Fig.  2  reverberation  smoothed  by  a  33  ms  rectangle 
from  a  91  m  SUS  charge  for  this  site  is  illustrated  in  the 
100Hz-2kHz  band.  This  data  was  collected  during  the 
SCARAB  98  trial  conducted  by  SACLANTCEN  (now 
NATO  Undersea  Research  Centre  (NURC)).  Clearly 
evident  is  the  dominant  scattering  from  the  region  of 
shoaling  bathymetry  and  thinning  sediment  cover  associated 
with  the  Ragusa  ridge.  Also  evident  are  strong  scattering 
events  from  certain  known  wrecks  to  the  South  and  East  of 
the  source  and  the  Campo  Vega  oil  production  platform  to 
the  NNW  (approximate  locations  are  indicated  by  the 
superimposed  black  symbols). 
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Figure  2.  Reverberation  data  collected  on  the  Malta  Plateau. 

We  are  interested  in  the  characteristics  of  clutter  in 
shallow  water  reverberation  datasets  such  as  illustrated  in 
Fig.  2.  Clearly  there  are  at  least  two  types  of  clutter  in  this 
data  set:  diffuse  clutter  characterized  by  increased 

reverberation  from  the  exposed  ridge,  and  discrete  clutter 
events  associated  with  anthropogenic  sources  such  as 
shipwrecks  and  oil  productions  facilities.  In  Fig.  3  we  show 
the  spectrogram  of  the  return  seen  on  the  bearing  of  1 1 6  true 
^  from  Fig.  2.  Here  one  can  observe  the  broadband  nature  of 
^  both  the  discrete  clutter  event  seen  at  15  s  and  the  extended 
!  return  from  the  Ragusa  ridge  from  19  to  24  s.  One  also 
!  observes  many  other  smaller  clutter  events,  most  energetic 
|  between  1  and  2  kHz.  Superimposed  on  these  features  is 
the  time-frequency  structure  of  coherent  round-trip 
propagation  to  the  scatterers  in  the  waveguide,  characterized 
by  a  clearly  visible  striation  pattern. 


Scarab  98  reverberation  116  true  (normalized  by  peak  at  1.33  s) 
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Figure  3.  Time-frequency  spectrogram  of  SUS  reverb  looking  towards 
shipwreck  and  Ragusa  ridge. 
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B.  Environmental  Description 

A  detailed  geoacoustic  and  scattering  description  of  the 
environment  between  Sicily  and  Malta  has  been  prepared 
based  on  measurements  of  propagation,  scattering  strength 
and  reverberation.  The  results  are  shown  in  Table  1.  Notice 
that  both  the  sediment  attenuation  and  the  scattering 
strength  grow  linearly  with  frequency.  Use  of  the  values  in 
Table  1  give  good  agreement  with  the  background  levels  of 
reverberation  on  the  sediments  of  the  Malta  Plateau,  away 
from  the  Ragusa  ridge  [7]. 


TABLE  I 

Background  Environmental  Properties  for  Malta  Plateau 


Property 

Value 

Sediment  Sound  Speed  Ratio 

1.05 

Sediment  Density 

1.80  gm/cmA3 

Sediment  Attenuation  (dB/A) 

0.70  f  (fin  kHz) 

Sediment  Scattering  Strength 
(Lommel-Seeliger) 

-47+14.8  f (fin kHz)  dB 

Ridge  Scattering  Strength 
(Lommel-Seeliger) 

-15  dB 

IV.  Simulation  Results 


Figure  4.  Comparison  of  data  and  prediction  101  true  337.5-362.5  Hz. 

Using  the  environmental  parameters  shown  in  Table  1 
and  the  system  parameters  of  beamwidth  and  source  depth, 
we  exercised  the  R-SNAP  model  to  predict  the  expected 
value  of  the  sediment  and  basement  reverberation  for 
various  beam  angles  over  various  bands.  The  results  for  the 
beam  azimuth  of  101  true  and  the  337.5-362.5  Hz  frequency 
band  are  shown  superimposed  on  the  data  in  Fig.  4.  Notice 
that  the  increase  in  bottom  reverberation  over  the  ridge 
structure  for  the  sediment  interface  scattering  mechanism, 


shown  in  green,  is  quite  slight,  while  the  scattering  from  the 
basement  is  essentially  non-existent  until  the  thinning  of  the 
sediment  cover  over  the  ridge  allows  it  to  become  a 
significant  mechanism  at  a  time  of  approximately  18  s.  The 
offset  between  the  data  and  the  simulation  is  due  to  the  use 
of  insufficiently  accurate  estimates  of  sediment  thickness. 

Using  the  new  time  series  simulation  capability  of  R- 
SNAP,  we  simulated  reverberation  time  series  for  the  337.5- 
362.5  and  the  950-1050  Hz  bands  over  the  clutter-inducing 
Ragusa  ridge.  The  data-model  comparison  for  the  337.5- 
362.7  Hz  band  is  shown  in  Fig.  5,  and  for  the  950-1050  Hz 
band  in  Fig.  6.  In  each  figure  the  data  is  shown  in  blue  and 
the  simulation  in  red,  with  20  dB  removed  so  as  to  facilitate 
comparison.  Agreement  in  the  lower  band  is  quite  good, 
including  up  and  onto  the  reverberation  enhancement  on  the 
ridge.  In  the  upper  band,  the  agreement  is  less  good  as  the 
sediment  thickness  causes  a  significant  reduction  in  the 
contribution  of  the  basement  scattering  contribution.  In  Figs. 
7  and  8  the  spectrograms  of  the  data  and  the  model  are 
shown  side  by  side  for  the  two  bands.  Note  the  time- 
frequency  striations  visible  in  the  data  are  replicated  in  the 
simulations.  Many  of  the  aspects  of  the  clutter  caused  by  the 
Ragusa  ridge  are  faithfully  replicated  in  the  simulations. 

V.  Conclusions 

The  R-SNAP  code  has  been  modified  to  simulate 
reverberation  time  series,  and  the  code  has  been  exercised  to 
model  clutter  caused  by  an  exposed  ridge  on  the  Malta 
Plateau.  Model-data  comparisons  at  350  and  1000  Hz  show 
favorable  agreement  between  the  reverberation  data  from 
the  ridge  and  the  simulations  using  average  environmental 
parameters  supplied  by  an  independent  analysis.  These 
results  illustrate  that  with  sufficient  environmental 
knowledge  certain  clutter  characteristics  of  shallow  water 
reverberation  associated  with  large  environmental  features 
may  be  predicted. 
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RSNAP  synthetic  vs  data  (Scarab  98,  950-1050  Hz  (36°22.87  N  14-45.21  E  101  °  true) 
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Figure  5.  Data  (blue)-model  (red)  comparison  in  the  337.5-362.5  Hz  band 


Figure  6.  Data  (blue)-model  (red)  comparison  in  the  950-1050  Hz  band 


RSNAP  simulation  332.5-367.5  Hz  (36-22.87  N  14-45.21  E  101  -  true) 
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Scarab  98  data  332.5-367.5  Hz  (36  -22.87  N  14-45.21  E  101  -  true) 
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Figure  7.  Spectrogram  of  350  Hz  simulation  (left)  and  data  (right) 


RSNAP  simulation  950-1050  Hz  (36-22.87  N  14-45.21  E  101  -  true) 
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Figure  8.  Spectrogram  of  1kHz  simulation  (left)  and  data  (right) 


